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Abstract 

The gravitational collapse of a magnetised medium is investigated by studying qualitatively the 



> 

IT) 

convergence of a timelike family of non-geodesic worldlines in the presence of a magnetic field. Focus- 
ing on the field's tension we illustrate how the winding of the magnetic forcelines due to the fluid's 
rotation assists the collapse, while shear-like distortions in the distribution of the field's gradients resist 
contraction. We also show that the relativistic coupling between magnetism and geometry, together 
with the tension properties of the field, lead to a magneto-curvature stress that opposes the collapse. 
This tension stress grows stronger with increasing curvature distortion, which means that it could 
potentially dominate over the gravitational pull of the matter. If this happens, a converging family of 
non-geodesic lines can be prevented from focusing without violating the standard energy conditions. 



PACS number(s): 04.20.-q, 04.40.-b 

1 Introduction 

Magnetic fields are very common in astrophysical environments and stellar magnetism is a long established 
and very active branch of astrophysics. Nevertheless, the study of magnetic, and of electromagnetic, 
fields in strong gravity environments is less developed. Most of the available studies address the possible 
gravitational effects on the Maxwell field and relatively few look into the implications of magnetic fields, 
in particular, for gravitational collapse itself. One of the most intriguing results so far has been obtained 
by Thorne in his analysis of Melvin's cylindrical magnetic universe [1] . There, by developing the concept 
of 'cylindrical energy', the author reached the conclusion that 'a strong magnetic field along the axis of 
symmetry may halt the cylindrical collapse of a finite cylinder before the singularity is reached' [2]. The 
possible support of the Maxwell field against the gravitational collapse of massive bounded systems was 
also studied in [3] . That work led to solutions of the Einstein-Maxwell system where the gravitational 
attraction is solely balanced by magnetic stresses. Studies of contracting charged dust have also suggested 
that the fluid may 'rebounce', thus preventing black-hole formation [4]. It has been argued, on the other 
hand, that a collapsing spherically symmetric charged dust will produce naked singularities due to shell- 
crossing [5]. Although these singularities are considered weak because the curvature invariants and 
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the tidal forces remain finite [6], their appearance could also signal that a nonzero Lorentz force and 
spherically symmetric charged collapse may not be physically compatible [7]. 

In this article we revisit the issue of the magnetic impact on gravitational collapse from an apparently 
entirely different viewpoint. We consider the non-spherical (but not necessarily cylindrical) collapse of a 
magnetised fluid, by studying the convergence of two neighbouring particle worldlines. As it turns out, we 
arrive at the same qualitative result as Thorne did, namely that the gravitational collapse of a magnetised 
fluid may stop before reaching the singularity. In our study, however, the reasons are seemingly unrelated 
to the cylindrical energy of Melvin's universe or to the charge density of the magnetised matter. We find 
instead that it is the intricate coupling between magnetism and geometry and the tension properties 
of the magnetic forcelines, namely their elasticity, which gives rise to ever increasing resisting stresses 
that may prevent the ultimate collapse from happening. Although these magneto-geometrical tension 
stresses are nothing more than the relativistic generalisation of a rather well known Newtonian effect, 
their existence remains largely unknown in the literature. Among the effects of the tension stresses we 
identify what is usually referred to as 'magnetic braking' and demonstrate how it assists the collapse of 
the magnetised fluid. We also find that shear-like distortions in the distribution of the field's gradients 
resist contraction. 

Studying the gravitational collapse of a magnetised fluid in full means solving the complete Einstein- 
Maxwell equations; a formidable task even for the most powerful numerical codes [8, 9]. Given that, 
all analytical studies so far have allowed for a considerable degree of mathematical simplification. Here 
we will consider the question of magnetised collapse without attempting to solve the Einstein-Maxwell 
equations. The main aim of this qualitative study is to shed light upon the role of the magnetic tension and 
draw attention to its potential implications. We assume a highly conductive fluid, similarly to the recent 
numerical work of [10], but no axial symmetry. This difference could be the reason the aforementioned 
numerical work did not encounter the resisting magneto-curvature effects discussed here. On the other 
hand, it could also be that these tension stresses dominate the dynamics of the collapse under special 
circumstances only. 

2 The worldlines of magnetised matter 

We will study the magnetic implications for gravitational collapse qualitatively by testing the convergence 
of the particle worldlines. In doing so we will be using the covariant approach to general relativity [11], and 
in many respects our discussion will resemble that of [12]. Throughout this article we assume conventional 
matter with positive gravitational mass and pressure, which means that the standard energy conditions 
are always fulfilled. High conductivity means that there is no electric field and that the magnetic field is 
'frozen in' with the matter. This is the well known MHD approximation [13]. 

When looking into the dynamics of gravitational collapse Raychaudhuri's formula is the key equation, 
as it covariantly describes the volume evolution of a self-gravitating fluid element [14]. Consider a 
congruence of timelike worldlines tangent to the 4- velocity field u a (with u a u a = —1). These are the 
worldlines of the fundamental observers and follow the motion of the fluid. The Raychaudhuri equation 
determines the proper-time evolution of © = V a u a , the scalar measuring the average contraction (or 
expansion) between two neighbouring worldlines [11]. In a magnetised environment we have [15] 

G = -\® 2 - \k (h + 3p + B 2 ) - 2a 2 + 2to 2 + B a ii a + ii a u a , (1) 

where k = 8irG, p and p are respectively the energy density and pressure of the fluid, B 2 = B a B a 
measures the energy density and the isotropic pressure of the magnetic field (B a ), a 2 and u) 2 are the 
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respective magnitudes of the shear and the vorticity associated with u a and ii a = u b VbU a is the 4- 
acceleration. The latter satisfies the momentum-density conservation law, which for a highly conductive 
perfect fluid takes the form [15] 

(H + p + \B 2 ) u a = -B a p - e abc B b cui\B c - U ab ii b , (2) 

where D a = h a b V a is the covariant derivative operator orthogonal to u a and Ii a b = —B^B^ describes 
the magnetic anisotropic pressure 1 . Note that we consider non-geodesic worldlines, since the motion of 
the particles is dictated by the combined Einstein-Maxwell field and not by gravity alone. Also, the fluid 
flow is not hypersurface orthogonal which explains the presence of the vorticity term in (1). 

The right-hand side of Eq. (1) determines the dynamics of the average volume evolution. Terms 
that are positive definite lead to expansion, while negative definite terms cause contraction. Thus, when 
the standard energy conditions are satisfied, all the right-hand side terms have a clear-cut role with the 
exception of B a ii a , which in principle can go either way. For the rest of this study, we will focus on B a ii a 
and examine its potential implications for the final fate of a collapsing magnetised fluid element. 

3 Lorentz force and magnetic tension 

The magnetic contribution to D a ti a comes from the Lorentz force. The latter is always normal to the 
direction of the field lines and emerges whenever the magnetic pattern is distorted from the condition of 
local equilibrium. Here, the Lorentz force is determined by the acceleration vector 

a a = -e abc B b cur\B c = -\B a B 2 + B b B b B a , (3) 

with a a B a = 0. The first gradient in the right-hand side of the above is due to the magnetic pressure and 
the second comes from the field's tension (e.g. see [13]). Insofar the tension stress is not balanced by the 
gradients of the field's pressure, a net force is exerted on the fluid particles. Note that e abc is the totally 
antisymmetric alternating tensor orthogonal to u a and therefore a a u a = 0. Using the 3-Ricci identity 
(see Eq. (6) below), the projected gradient of a a decomposes as 

B b a a = — \B b B a B 2 + B b B c B c B a + B c B c B b B a + K acbd B c B d - 2u bc B {a) B c . (4) 

Here TZ a bcd is the Riemann tensor of the observer's instantaneous 3-dimensional rest space and u> ab is the 
vorticity tensor associated with the fluid flow. The last four terms on the right-hand side of the above 
convey the magnetic tension effects. In particular, the magneto-curvature stress in Eq. (4) reflects the 
special status of vectors, as opposed to that of scalars, in general relativity. This special status stems 
from the geometrical nature of the theory and it is manifested in the Ricci identity 

2V [a V b] B c = R dcba B d , (5) 

applied here to the magnetic field vector, with R a bcd being the spacetime Riemann tensor. When projected 
orthogonal to u a , the above expression leads to what is commonly referred to as the 3-Ricci identity [11] 

2B [a B b] B c = n dcba B d -2uo ab h c d B d . (6) 

1 Angled brackets indicate the symmetric, trace-free part of projected second-rank tensors (e.g. B^ a B b ) = B a B b — 
(B 2 /3)hab, where h ab = g ab + u a u b with h ab u b = 0). 
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The Ricci identities argue for a direct interaction between vector sources and spacetime curvature, which 
adds to the standard interlay between matter and geometry as we know it from the Einstein field equa- 
tions. In the magnetic case this direct coupling also brings into play the tension properties of the latter, 
namely the elasticity of the magnetic force lines, and couples it in an intricate way with the geometry of 
the space. This unique feature, the magnetic tension, manifests the well known fact that the field lines 
'do not like to bend' and react to any attempt that distorts them. Indeed, within the Newtonian theory 
the magnetic tension is known to trigger restoring stresses which depend on the strength of the field and 
on the deformation of the magnetic force lines (measured by their curvature radius) [13]. What the gen- 
eral relativistic expression (4) shows, is that deviations from Euclidean geometry will lead to analogous 
tension stresses, which are also proportional to the magnetic strength and to the amount of curvature 
distortion. This time, however, it is the curvature of the space itself that causes part of the magnetic 
deformation. In other words, the magnetic field 'feels' the curvature of the space in a way that is dictated 
by the Ricci identities and this is demonstrated by the magneto-curvature stress in the right-hand side 
of Eq. (4). The effects of this tension stress are generally counter- intuitive because of the nature of 
the magnetic property itself and its subtle coupling with the geometry of the space. The latter means 
that even weak magnetic fields can lead to a strong overall effect under favourable circumstances. The 
potentially pivotal implications of magnetism, and of the magnetic tension in particular, for cosmology 
were originally discussed in [16]. These results have been confirmed and extended in [17], although an 
oversight prevented the authors of the latter paper from recognising the key role of the magnetic tension. 
Here, we will consider the implications of the same tension stresses for the gravitational collapse of a 
magnetised fluid. 

Suppose that both the fluid and the magnetic field have a nearly homogeneous energy density distri- 
bution. In practise this means that any inhomogeneities that might be present are relatively small and 
that we may ignore spatial gradients in the energy density of the two sources. For the magnetic field 
our assumption implies that the Lorentz force term in the right-hand side of (2) is dominated by the 
tension stresses. This is exactly what we need, since our aim is to investigate the effect of these particular 
magnetic stresses. Then, Eq. (2) gives u a B a = and subsequently it reduces to 

+ p + B 2 ) ii a = -e abc B b cm\B c , (7) 

Taking the projected divergence of the above, using the trace of (4), with T) a B a = due to the absence 
of electric fields, and substituting into Eq. (1) we arrive at 

9 + ±9 2 = -\k (fi + 3p + B 2 ) -2 (a 2 - £ 2 ) + 2 (to 2 - W 2 ) + e- 1 K ab B a B b + u a ii a , (8) 

where e = n+p + B 2 , S 2 = (D {a B b) ) 2 /2e, W 2 = (D [a B b] ) 2 /2e and K ab = TZ c acb is the 3-Ricci tensor given 
by the Gauss-Codacci equation [11, 18] 

K ab = h a c h b d R cd + R acbd u c u d - k c c k ab + k ac k c b , (9) 

where k ab = B b u a is the extrinsic curvature tensor. The second term in the right-hand side of the above 
ensures that, in addition to the usual matter fields, the Weyl curvature (i.e. tidal forces and gravitational 
waves) also contributes to the geometry of the 3-space. Also note that in deriving Eq. (8) we have 
assumed magnetic flux conservation, which guarantees that the rotation term in the trace of (4) vanishes. 
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4 The magnetic tension stresses 



We will first turn our attention to the quantities E 2 and W 2 in the right-hand side of (8), which are 
the magnitudes of E a b = D^B^/^/e and W a b = Ti^B^/^/e respectively. The former of these tensors 
describes distortions in the distribution of the field gradients and resists contraction, while the latter is 
the magnetic twist tensor and assists the collapse. Note that, although E a & and W a & have the shear and 
the vorticity as their respective kinematical analogues, their effect on is exactly the opposite of the 
one normally associated with the shear and the vorticity proper (see also [19]). This counter- intuitive 
behaviour reflects the fact that both E and W carry the tension properties of the field and manifests the 
tendency of the magnetic force lines to remain straight. Thus, while a nonzero shear assists the collapse, 
the corresponding magnetic stress tries to balance this effect out. Consider also the tensor W a b, which 
has magnitude W 2 = W a bW ab /2 = (curLB a ) 2 /4e and is triggered by the winding of the magnetic field 
lines around a rotating fluid element. Following (8), a nonzero W will always reduce the gravitational 
effect of kinematic vorticity. This manifestation of the field's tension is also known as 'magnetic braking' 
and can accelerate the gravitational collapse of a rotating star [10]. 

For our purposes the key quantity is the magneto-curvature tension stress lZ a i } B a B b in the right- 
hand side of Eq. (8), which measures the curvature of the 3-space along the direction of the magnetic 
force lines. Starting from (9) one can show that, when only the magnetic field is present, lZ a \ ) B a B b = 
[20]. Therefore, despite the magnetic energy input, the curvature of the 3-space in the direction of 
the field lines is zero; a result independent of the magnetic strength. Technically speaking, it is the 
negative pressure of the field along its own direction which cancels out the positive contribution of the 
magnetic energy density. More intuitively, it is the inherit tendency of the magnetic force lines to remain 
'straight' which is responsible for the aforementioned null result. Clearly, in the presence of other sources 
lZ a bB a B b / 0. Then, the magneto-curvature effect on is rather unexpected. For 1Z a i ) B a B b < the 
magneto-geometrical stress in (8) brings the particle worldlines closer, but pushes them apart when the 
field lines are 'positively curved', that is for lZ a t,B a B b > (see [16] and also [17]). This is against 
the common perception, which always associates positive curvature with gravitational contraction. As 
with the magnetic shear and the magnetic vorticity stresses discussed earlier, the reason for the counter- 
intuitive behaviour of the magneto-curvature term is the tension properties of the field lines. 

Let us assume, mainly for illustration purposes, that the effect of E and W is cancelled out by that 
of their kinematic counterparts. In other words, we will consider the case where uj 2 + E 2 = W 2 + a 2 . 
Although it may not appear so initially, we will later show that this assumption is much less restrictive 
that it looks. For the moment we note that, in addition to counteracting each other, the pairs uo 2 , W 2 and 
a 2 , E 2 are of the same nature and 'perturbative order' (i.e. quadratic in D a Ub and D a Bb). The former of 
these properties supports the assumption that the aforementioned opposing pairs are very likely to balance 
each other out. The latter ensures that these terms, unlike H a \ ) B a B b for example, become appreciable 
only in highly inhomogeneous configurations. Under such conditions, the Raychaudhuri equation reads 

+ |0 2 = -R ab u a u b + c 2 U ah n a n b + u a u a , (10) 

where n a = Ba/^/B 2 , c 2 = B 2 /e is the Alfven speed and R a b is the Ricci tensor of the spacetime with 
R ab u a u b = k(h + 3p + B 2 )/2 > 0. The latter means that the strong energy condition is satisfied. 

It should be noted that, in addition to the aforementioned arguments regarding the relative gravi- 
tational input of the shear, the vorticity and of their magnetic counterparts, bypassing these terms also 
helps to isolate the two resisting stresses in Eq. (8), namely c 2 1Z a b l o. a n b and u a u a - The latter has been 
analysed within a collapsing perturbed Tolman-Bondi spacetime in [7], where it was found to grow faster 
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than all the other terms in the right-hand side of the Raychaudhuri equation. Here we will focus on the 
magneto-curvature tension stress instead. In any physically realistic scenario of stellar collapse this term 
is always positive and therefore it always resists against further contraction. Moreover, the strength of the 
gravito-magnetic stress is essentially proportional to the amount of the curvature distortion. This feature 
distinguishes the tension stress from the rest and makes it a very promising candidate for outbalancing 
the gravitational pull of the matter. 

5 Gravitational pull vs gravito-magnetic tension 

The assumption of a spatially homogeneous energy density distribution for the sources means that D a p = 
0, which in turn guarantees that the 4-acceleration vector ii a depends entirely on the magnetic field (see 
Eq. (7)). Therefore, when the field is absent all the positive definite terms in the right-hand side of (10) 
vanish and an initially converging congruence (i.e. one with ©o < 0) will focus (i.e. — > -co) within a 
finite amount of time, unless the energy conditions are violated. This is a fundamental and well known 
result about gravitational collapse (e.g. see [18]). In the magnetic presence, however, the two positive 
definite terms in the right-hand side of (10) will resist the collapse. Thus, ignoring the supporting effect of 
u a u a , we argue that the magneto-curvature effects can prevent a converging congruence of non-geodesic 
worldlines from focusing, without violating the standard energy conditions, if 

clK ab n a n b > R ab u a u b . (11) 

Obviously, the above also holds when the earlier imposed constraint to 2 + E 2 = W 2 + a 2 is replaced by 
uj 2 + £ 2 > W 2 + a 2 . It still holds when uj 2 + S 2 < W 2 + a 2 , provided that u a u a /2 > W 2 + a 2 - to 2 - S 2 . 
So, condition (11) holds for a variety of combinations between uj 2 , a 2 , W 2 and S 2 . This means that 
the earlier restriction placed on these quantities is not essential for the validity of our argument. The 
same can also be said about the fluid inhomogeneities, given that pressure gradients resist the collapse 
(through the last term in the right-hand side of (10)). Our approximations, however, have helped to 
isolate and demonstrate the role of the curvature terms in Eq. (8), which should decide the ultimate fate 
of the collapse. Indeed, the main reason for focusing on the two geometrical quantities in the right-hand 
side of (8) is that, as the collapse proceeds, we expect the curvature to dominate. Then, the fate of the 
collapsing magnetised fluid should be decided by the balance between the two quantities in Eq. (11) and 
the curvature terms in (8) have not been unduly favoured at the expense of the rest. If the gravitational 
pull of the matter, which is the driving force behind the collapse, is outbalanced, there is a realistic 
possibility of avoiding worldline focusing. The magneto-curvature tension stresses open this possibility. 
This is so because, while only the usual matter fields contribute to R a bU a u b , the tension stress TZ a i,n a n b 
has additional contributions from other sources (see Eq. (9)). Probably the most important among these 
extra sources is the Weyl curvature. This monitors the long-range gravitational field and contributes to 
1Z a b directly via the electric Weyl component [11]. The latter describes the tidal forces which are expected 
to increase dramatically as the collapse proceeds. These additional contributions to the 3-Ricci curvature 
mean that, in principle, the magneto-curvature tension stresses can outbalance the gravitational pull of 
the matter. Although the final outcome depends on the precise characteristics of the collapse, our analysis 
points towards the conjectural but very intriguing possibility that if the magnetic line deformation due 
spatial curvature distortions is strong enough, the resulting tension stresses may just be able to avert the 
formation of caustics and eventually the ultimate collapse of the magnetised fluid. 
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6 Discussion 



The immediate implication of the above study is that violating the standard energy conditions to prevent 
an initially converging congruence from focusing may not be always necessary when magnetic fields are 
present. In particular, nonspherical magnetised collapse (the magnetic presence will inevitably distort 
spherical symmetry to a larger or lesser degree) may not end up in a 8 -> -oo singularity because of the 
field's tension. Whether our analysis applies to physical situations, such as the gravitational implosion of 
a massive star, depends on what the exact properties of a realistic collapse are. It might be, for example, 
that stellar magnetic fields do not survive to the later stages of the collapse, or that the MHD limit is 
no longer a good approximation. We have no reasons to believe that either of these possibilities may be 
true however. In fact, the opposite appears more likely given the presence of very strong magnetic fields 
in compact stellar objects like neutron stars. It might also be that the gravitational pull of the matter 
always prevails and condition (11) is never satisfied. Very recent numerical simulations of magnetised 
hypermassive neutron-star collapse indicate that, at least when axial symmetry holds, this seems to be 
the case [10]. If condition (11) is met, however, the magnetic tension could stop the convergence of the 
fluid worldlines and the same could also happen to the contraction of a magnetised star. Similarly, if 
one evolves a magnetised universe backward in time, they may find a highly curved state that could 
potentially re-expand into the past. It should be noted, however, that avoiding worldline convergence 
only means avoiding a singularity in the congruence and does not guarantee a singularity-free spacetime; 
at least under the current consensus of what a singularity is [18]. For example, the spacetime can still 
be geodesically incomplete. Nevertheless, any such singularities should be of limited influence if most of 
the matter can successfully avoid them. 

Stresses that support against gravitational collapse are not exclusively particular to a magnetic pres- 
ence. It is well known that whenever the worldlines are not hypersurface orthogonal or geodesies, sup- 
porting stresses always appear due to rotation or pressure gradients. Vorticity, for example, has been 
considered in the past as a possible way of preventing the ultimate collapse [21]. Therefore, it is not so 
much the existence of magnetic related stresses that resist gravitational collapse, as the nature of the 
stresses themselves. In this respect, our key result is that, when magnetic fields are involved, one of 
the supporting stresses depends (in fact it is proportional to) on the distortion of the curvature itself. 
The presence of this stress is a direct and inevitable consequence of the vector nature of magnetic fields 
and of the geometrical nature of general relativity, while their counter-intuitive effects result from the 
tension properties of the field. It is the elasticity of the magnetic force lines, their inherit tendency to 
remain 'straight', which manifests itself as a reaction to curvature distortions that is proportional to the 
distortion itself. In a sense, it appears as though the elastic properties of the field have been injected in 
to the fabric of the space. 

Finally, when considering the unconventional magnetic behaviour described so far, the reader should 
keep in mind that magnetic fields are rather unusual sources themselves. The Maxwell field is the only 
vector source that we know that exists in the universe today. Within the geometrical framework of general 
relativity, the status of vector fields is different from that of the ordinary matter. This is so because, in 
addition to the Einstein field equations, vector sources interact with the spacetime geometry through the 
Ricci identities as well. This purely geometrical coupling, which has been at the centre of our discussion, 
is already known to trigger some other rather nontrivial effects. The best known example is probably 
the 'scattering' of electromagnetic radiation by the gravitational field, which leads to what is commonly 
referred to in the literature as 'wave tales' [22]. Here, we have outlined the potential implications of 
effectively the same Einstein-Maxwell coupling for gravitational collapse. 
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